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Lately, it has been suggested that strange (quark) stars can be responsible for glitches and other 
observational features of pulsars. Some discussions on whether quark stars, if really exist, are bare 
or crusted are also a source of controversy in the recent literature. In the present work we use the 
Nambu-Jona-Lasinio model, known to incorporate chiral symmetry, necessarily present in the QCD 
formalism, in order to describe quark star properties. We compare our results for the stars and 
the features of the model with the much simpler model normally used in the description of strange 
stars, namely the MIT bag model. We also investigate the differences in the stellar properties which 
arise due to the presence of the crust. We show that the NJL model produces results which are 
somewhat different as compared with the MIT model. 
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, I. INTRODUCTION AND FORMALISM 

r> , 

Neutron stars are the remnants of supernova explosions with masses 1-2M Q , radii ~ 10 km, and a 
temperature of the order of 10 11 K at birth, cooling within a few days to about 10 10 K by emitting 
neutrinos. In a conventional model for a neutron star, the star is composed of hadrons, predominantly 
• degenerate neutrons with an admixture of protons, and degenerate electrons. In the stellar modeling, 

Oh' the structure of the star depends on the assumed equation of state (EOS), which contains a number of 

uncertainties associated with uncertainties in the strong force under conditions appropriate to a neutron 
£^ ■ star. An important uncertainty concerns the true ground star of matter. In conventional models, hadrons 

are assumed to be the true ground state of the strong interaction. However, it has been argued 0,0, Q 
that 'strange matter' is the true ground state of all matter. If this is correct, the interior of neutron stars 
should be composed predominantly of u, d, s quarks, plus leptons to ensure charge neutrality. This led to 
the suggestion that there may be no neutron stars, and that all neutron-like stars are in fact strange stars 
PJ. More generally, the composition of neutron stars remains a source of speculation, with some of the 
possibilities being the presence of hyperons a mixed phase of hyperons and quarks 0, IE 01 El El > 

a phase of deconfined quarks or pion and kaon condensates [l^. Models in which the interior is assumed 
to be composed of strange matter are often called 'strange' stars. However, because, as we show, the 
strangeness content depends on the model used to describe the quark matter, we prefer to describe any 
model in which the interior involves deconvolved quarks (not bound in hyperons) as 'quark' stars. In the 
stellar modeling, the structure of the star depends on the assumed EOS, which is different in each of these 
cases, and which depends on the nature of the strong interaction. Apart from the differences in the EOS, 
an important distinction between quark stars and conventional neutron stars is that the quark stars are 
self-bound by the strong interaction, whereas neutron stars are bound by gravity. This allows a quark star 
to rotate faster than would be possible for a neutron star. Further evidence in favor of quark stars is that 
some stars do indeed seem to rotate faster than what would be expected for a neutron star 0, EL El . 

A quark star must have a thin layer on its surface dominated by the electrons which are necessary to 
enforce charge neutrality. This layer could suspend a hadronic crust, which would not be in contact with 
the stellar core It has been suggested that the presence of a crust provides a natural explanation 

for pulsar glitches 23, which are sudden changes in the rotation period of the pulsar. Other authors have 
argued that quark stars should be 'bare' [Til Il7j . in the sense that any such crust would either not form 
or would be destroyed during the supernova explosion. Recently, the characteristics of the radiation from 
hot, bare strange stars have been identified |18| . 

Possible candidates for quark stars include compact star with gravitational masses around < IMq 
and radii of the order of ~ 6km: PSR J1645-0317 (PSR B1642-03), PSR J1830-1059 (PSR B1828-11) 
EH, RX J185635-3754 HJ, Her X-l [U, 4U 1728-34 gj, X-ray bursters GRO J1744-28 || and SAX 
J1808.4-3658 El. 
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There are several different models that describe quark matter. In |25l we investigated the properties 
of strange stars using two different models: the MIT bag model |26|. which is widely favored in the 
literature on strange stars, and the color-favor-locked phase (CFL) model [27|], which allows the quarks 
near the Fermi surface to form Cooper pairs which condense and break the color gauge symmetry psj . 
At sufficiently high density the favored phase is called CFL, in which quarks of all three colors and all 
three flavors are allowed to pair. In p5j we verified that quark stars have different maximum gravitational 
and baryonic masses, compared with neutron stars, and that the maximum masses are not significantly 
affected by the presence or absence of a crust. Moreover, with a simple prescription for the Kepler 
frequency, which determines the spin rate of a neutron- like star, we obtained much higher values for the 
quark stars than for the neutron stars. Finally, in [25j | we used the two possible mass-to-radius ratio 
constraints available in the literature [29ll30| to distinguish between the assumed EOSs. Concerning one 
of the constraints, specifically the interpretation of absorption features as atomic transition lines in |30j| 
is controversial: an alternative interpretation |3lL l32| is that the aborption features are cyclotron lines, 
which imply no obvious constraint on the EOS. However, if the accepted interpretation is correct, then 
only the quark star EOS is compatible with the constraint given in |29| . thereby excluding most of the 
neutron star EOSs. The conclusion that all neutron stars must be quark stars has major implications, 
but it is necessarily valid only for either of the two EOSs investigated in [2^| . 

In the present paper we focus on a different model for the strange matter, the Nambu-Jona-Lasinio 
(NJL) model |3^|. Our aim is to use this model to obtain the stellar properties and compare the results 
with the MIT bag model 26]. It is important to distinguish between the EOS during the short time period 
when neutrinos are still trapped in the star, and the EOS after the neutrons escape. The maximum ent ropy 
per baryon (S) reached in the core of a new born star is about 2 (in units of Boltzmann's constant) |3 I| . 
We then perform our calculations for S = 0(T = 0), 1 and 2 since the entropy and not the temperature 
should be constant throughout the star |35|. 

In section II the formalism used is presented; in section III we give the results and make the relevant 
discussions and in section IV the conclusions are drawn. 



II. QUARK MATTER MODELS 

In this section we summarize the main formulae for both models used in this paper. 



A. The Nambu-Jona-Lasinio Model 

We choose the NJL model [3^, H3, IH, to describe the quark phase. The SU(3) version of the model 
includes most of symmetries of QCD, including chiral symmetry, and its breaking, which is essential in 
treating the lightest hadrons. The NJL model also includes a scalar-pseudoscalar interaction and the 't 
Hooft six fermion interaction that models the axial U(1)a symmetry breaking. The NJL model assumes 
deconfined point-like quarks, and is not remornalizable, requiring regularization through a cutoff in three- 
momentum space. 

The NJL model is defined by the Lagrangian density 
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L = q(i-fd M - m)q+gs ^[(gA^) 2 + ( qi 75 A a q f ] 

+ 9D {det [qt (1 + 75) qj] + det [q t (1 - 75) qj] }, (1) 

where q = (u, d, s) are the quark fields and A a ( < a < 8 ) are the U(3) flavor matrices. The model 
parameters are: m — diag (m u , , m s ) , the current quark mass matrix (m^ = m u ), the coupling 
constants gs and go and the cutoff in three-momentum space, A. The NJL model is valid only for quark 
momenta smaller than the cut-off A. 

The set of parameters is chosen in order to fit the values in vacuum for the pion mass, the pion 
decay constant, the kaon mass and the quark condensates. We consider the set of parameters [39L l40j: 
A = 631.4 MeV, g s A 2 = 1.824, g D A 5 = -9.4, m u = m d = 5.6 MeV and m a = 135.6 MeV which are 
fitted to the following properties: m, = 139 MeV, /„. = 93.0 MeV, m K = 495.7 MeV, f K = 98.9 MeV, 
(uu) = {dd) = -(246.7 MeV) 3 and (ss) = -(266.9 MeV) 3 . 

The thermodynamical potential density is given hy ft = £ — T S — SiMi-^j — ^0, where the energy 
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density is 

£ = -2N C ^J j—^ (m--n i+ )9(A -p) 

- 2g s (^q t ) 2 ~2g D (uu){dd)(ss) -£ Q (2) 

i=u,d,s 

and the entropy density is 

S = -2N C J] J -^9(A 2 -p 2 ){[n l+ \n(n l+ ) + (l-n l+ )\n(l - + [n i+ - ra^]} . (3) 

z— u,d,s 

In the above expressions AT C = 3, T is the temperature, /Uj (JV<) is the chemical potential (number) of 
particles of type i, and £o an d are included in order to ensure £ = O = in the vacuum. This 
requirement fixes the density independent part of the EOS. The ground state of the system is described 
by the density matrix jj^j given by / = diag (f u , fd, f s ) with 

f i = \[I{n i -+n i+ ) l0Mi + a ' P ( ni --n i+ )]6(A 2 -p 2 ), (4) 

where / is the identity matrix, n\ are the Fermi distribution functions of the negative (positive) energy 
states, n\ = [l + exp(=p(/3 (Ei±fXi)))]~ x , i — u,d, s. In the last equation /3 = 1/T, M, is the constituent 
quark mass, Ei — (p 2 + M 2 ) 1 ! 2 . 

The quark condensates and the quark densities are defined, for each of the flavors i — u,d, s, respectively, 
as: 

(q iqi ) = -2N c J -0^^L( ni _ -n l+ )8(A 2 -p 2 ), (5) 



Pi = (ft f ft) = 2iV c j (ni_ + n l+ - 1) 6{A 2 - p 2 ). (6) 

Minimizing the thermodynamical potential Q with respect to the constituent quark masses Mi leads to 
three gap equations for the masses Mj 

Mi = m l -4gs(ftft) - 2 g D (qj qj)(q~k Qk) , (7) 

with cyclic permutations of i, j, k. 

We introduce an effective dynamical Bag pressure |4l) . 

B = 2N C £ J ^(Vp T +M~ 1 -Vp T +^)0{A 2 -p 2 ) 

i=u,d,s 

- 2 9s Y (ft*) 2 - 4g D (uu)(dd)(s s). 

i=u,d,s 

In terms of this quantity the energy density J2) takes the form 

£ = 2iV c ^ J AP^/ p T^ M -.^ n . + _ n ._ + ^ d (A 2 ~p 2 ) +B eff , B eff =B -B, (8) 

j— u,d,s 

where Bq = B Pu=Pd=Pv= Q. Writing the energy density in terms of B e ff allows us to identify this contri- 
bution as a Bag pressure and establish a relation with the MIT Bag model (3^, El discussed in the next 
section. 
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B. The MIT Bag Model 

The MIT Bag model [2^j has been extensively used to describe quark matter. In its simplest form, the 
quarks are considered to be free inside a Bag and the thermodynamic properties are derived from the 
Fermi gas model. The energy density, the pressure and the quark q density, respectively, are given by 



q=u,d,s 

p = i E / d p 1 9 p4 = (/«+ + /«-) - Ba ^ ( 10 ) 

? J Jp 2 +m 2 



(11) 



where 3 stands for the number of colors, 2 for the spin degeneracy, m q for the quark masses, Bag represents 
the bag pressure and the distribution functions for the quarks and anti-quarks are the Fermi distributions 



/ 9 ± = l/(l + exp[(e T ^)A1), (12) 

with p q (—fi q ) being the chemical potential for quarks (anti-quarks) of type q and e = (p 2 + m 2 ) 1 / 2 . 
These equations apply at nonzero temperatures. For T = 0, there are no antiparticles, and the particle 
distribution functions become the usual step functions. 

We use m u — rrid — 5.5 MeV, m s — 150.0 MeV and Bag = (180 MeV) 4 . If m u , and m s are chosen as 
in the NJL model, the behavior of the properties of interest are not altered, since they are more dependent 
on the Bag pressure than on small differences in the quark masses. 



C. Quark matter in beta equilibrium 

In a star with quark matter we must impose both beta equilibrium and charge neutrality . In what 
follows two scenarios are investigated, an early stage when there are trapped neutrinos in the interior 
of the star and a later stage, after the neutrinos escape (deleptonization). We first consider the later 
stage when entropy is maximum and neutrinos diffuse out. The neutrino chemical potential is then zero. 
For /3-equilibrium matter we must add the contribution of the leptons as free Fermi gases (electrons and 
muons) to the energy and pressure. The relations between the chemical potentials of the different particles 
are given by 

Ha = /id = PL U + He, He = /V (13) 

For charge neutrality we must impose 

1 . 

Pe + Pfi = g( 2 P« - PA - Ps)- 

For the electron and muon densities we have 

/TO 
j^(fl+-fl-), l = e,p, (14) 

where the distribution functions for the leptons are given in eq. (|12|l by substituting q by I, with p,i the- 
chemical potential for leptons of type I. At T = 0, eq. (fHEj) becomes pi = kpJSn 2 . The pressure for the 
leptons is 

p * = Ae/ K d l 2 (Si + +Si-). (15) 

In earlier stage, when the neutrinos are still trapped in the interior of the star, eq. (|13|l is replaced by 

Ps = Pd = Pu + Pe - Pue, (16) 

the lepton contribution is set to be Yl = Y e + Y ve = 0.4 |3f|. No muons appear in this case. 
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III. RESULTS 

In order to describe the crust of the quark stars we next use the well known EOS calculated in f° r 
very low densities. 

In figure ^ the strange quark mass as a function of the density is shown for the NJL models and the 
6 cases discussed in the present work, i.e., 5 — 0,1,2 without neutrinos and with trapped neutrinos. 
These data have a substantial influence on the results we obtain with the NJL model. Chiral symmetry 
restoration occurs first for the u and d quarks and only at higher densities for the s quarks. This is clearly 
seen in figure^ the mass of the s-quark starts to decrease only for p/ po > 4 (6) for neutrino free (trapped) 
matter. 

In figure [21 we show the EOS obtained respectively for the MIT bag and for the NJL models again 
for the 6 different cases mentioned above. As shown in the upper panel in figure |3 the results are very 
similar and it is difficult to distinghish between them; this may be attributed to the familiar result that 
for an extreme relativistic gas the pressure is one third of the energy density. The NJL model presents 
different behaviors at high densities. The fact that chiral symmetry is restored at different densities for 
quarks with different masses explains the softer behavior of the EOS for e > 4fm~ 4 . At these densities 
the onset of strangeness occurs, with the onset occurring at larger densities when neutrino trapping is 
included (see Fig.^J. The EOS with trapped neutrinos is harder and the S = 0, 1, 2 cases are similar, but 
distinghishable since higher temperatures correspond to harder EOSs. This behavior was already seen in 
H3, although in a different context. 

In figure |3| we show the temperature range described by both models for the cases 5=1 and S = 2. 
For the MIT model, the temperature for each different EOS increases steadily, reaching about 35 MeV for 
5 = 2 and ~ 17 MeV for S = 1 at p = 12po- The NJL model presents a different behavior: the temperature 
first oscillates before increasing. The near-plateau around 4-6 p/po coincides with the onset of strangeness 
opening of a new degree of freedom. The maximum temperature for 5 = 2 without neutrinos is higher 
than the one obtained with the MIT model for the same energy density. This difference is due to the 
fact that within the MIT model the fractions of the different types of quarks are similar, unlike the NJL 
which has a smaller fraction of s-quarks, which gives rise to a higher entropy. If we fix the entropy the 
temperature must be lower within the MIT model. However, from Table I we can see that within the NJL 
the central energy densities of compact stars are much lower than the corresponding ones for the MIT 
model. The temperatures are not shown in Table I because they vary though the star; the temperatures 
in the interior of the quark stars are then similar for the NJL and the MIT models, despite the significant 
difference in density. Including trapped neutrinos lowers the temperature of the star. This is due to the 
presence of a fixed fraction of leptons increasing the degrees of freedom, and then to keep the entropy per 
particle and the thermodynamic potential fixed, a lower temperature is required. If neutrino trapping is 
not imposed the fraction of leptons in the star is very small. 

In figure^we plot the strangeness content in each situation described in the text. As already referred, 
within the NJL the onset of strangeness occurs for p/po ~ 4 depending on the temperature and the neutrino 
content. A different situation applies for the MIT bag: strangeness is present from lower densities and at 
p/ Po ~ 4 the strangeness fraction is almost at its maximum value in the case without neutrinos, ~ 0.33 . 
The presence of leptons in the case of trapped neutrinos with a fixed lepton fraction lowers the strangeness 
content due to the electric charge conservation restrictions. 

In figure we display, for both the MIT and the NJL models, the electron neutrino fraction present in 
the interior of the stars when trapped neutrinos are included. With the MIT model the neutrino content 
is practically the same for the three entropies considered. It increases a little for lower densities and then 
reaches its maximum value, for the same reasons discussed above, i.e., neutrinos and electrons together 
are fixed to a certain fraction and electrons compensate for charge neutrality in a system where the three 
quarks have equal number densities. The NJL model shows a slight decrease in the neutrino fraction until 
the onset of strangeness. This effect is more pronounced for p < 3po when chiral symmetry restoration 
has still not occured for the u and d quarks. A slight increase in the neutrino fraction occurs with the 
increase of the entropy, mainly at low densities. The neutrino fraction never reaches 14% within the NJL 
picture, but it amounts to more than 15% within the MIT model. We see a similar neutrino fraction for 
these models in [n| . The difference in neutrino content is again due to the larger fraction of strange 
quarks in the MIT description and therefore a smaller number density of electrons required by the electric 
charge neutrality. A smaller fraction of electrons implies a larger fraction of neutrinos for a fixed lepton 
fraction. Nevertheless, if the CFL model were used for the quark matter with trapping, a much higher 
neutrino fraction, of the order of 35%, would be possible [Toj . 

Given the EOS, the next step is to solve the Tolman-Oppenheimer-Volkoff equations 01 . Table I lists 
the main properties of the stars, including their maximum gravitational and baryonic masses, their radii 
and central energy densities, for all the EOSs considered in the present work. As expected, quark stars 
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have smaller maximum masses with smaller radii 0, ITol l25j and higher central densities. 

Let us first analyze the results obtained with the MIT model and then outline the main differences that 
appear with the NJL model. For quark stars either without or with neutrinos, the maximum gravitational 
masses and radii are always around 1.23-1.24 Mq and 6.7-7.1, respectively, and the energy density is 
around 14.5 fm~ 4 . The inclusion of trapped neutrinos and a crust, have practically no effect on the 
numerical results. 

The NJL model leads to significantly different results. The maximum masses and radii are systematically 
larger and, as a consequence, the central energy densities are lower than for the MIT bag model. In fact, 
the results do not depend strongly on the entropy, as already pointed out in 25]. In the NJL, for a bare 
star without neutrinos, the maximum masses and radii decrease and the central energy density increases 
as the entropy increases. For a star with a crust, the gravitational and baryonic masses decrease but the 
radii oscillate as the entropy increases. 

In Fig. H3 we show for S = 0, and for both models the mass-radius plot of the family of stars obtained 
for the neutrino free and the neutrino trapped situations. It is clear that the NJL predicts the possibility 
of existence of more massive and with larger radius stars. For a fixed radius smaller than the larger 
accepted by the MIT model the NJL gives stars with much smaller gravitation masses. Conclusions for 
finite entropies are the same. 

Another difference is that when trapped neutrinos are included, this generally increases the maximum 
gravitational mass, with the increase being much larger for the NJL than for the MIT model. 

In a previous work [25[ we see that the MIT results are sensitive to the Bag parameter. With a lower 
value of this parameter, results similar to the ones obtained with the NJL model are found. The same is 
possible when the CFL model is used. In [25j we fixed the gap parameter of the CFL model as 100 MeV 
and showed that the properties of the star depend on the Bag parameter, but the properties are also 
known to depend on the gap parameter [Sj. 

It is not clear that quark stars really exist, and if the do exist in principle, how they would be formed. 
They form directly in a supernova explosion, or as a consequence of the hadron-quark deconfinement 
phase transition in stellar compact stars |45j . In the former case there is neutrino trappin g du ring the 
first seconds of the stars life and the entropy is approximatly constant throughout the star [35j. In this 
case one concludes from our results that there can be no late blackhole formation. This event is only 
possible if the baryonic mass of a stable star at finite temperature is larger than the largest baryonic mass 
of a stable cold star [34j]. From Table I one can see that the maximum baryonic masses at fixed entropies 
are always lower than the corresponding values at zero temperature. This kind of effect could occur within 
hybrid stars [Hill. 

From figure and Table I we can also see that the inclusion of a crust has a similar effect for both 
the MIT and NJL models: the radius of the maximum mass star increases ~ 0.3 km. It has been argued 
that at finite temperature the star tends to be bare due to the reduction of the electrostatic potential 
of the electron |4y|. In particular, at T = 30 MeV a star has essentially no crust. Taking into account 
simulations with neutron stars, we should consider the entropy rather than the temperature to be constant 
throughout the star |3{|. In this case the temperature at the surface is lower: we get ~ 8 and 16 MeV 
respectively for S — 1 and S = 2, respectively. 

An important difference between the NJL stars and MIT stars is that the electron chemical potential 
is very different. In Fig. we show, as a function of the density, the electron chemical potential in both 
models for the different values of entropy. Entropy has almost no effect on the results but there is an 
important difference between the MIT and the NJL cases: fi e < 20 MeV for the MIT model and ^ e as 
high as 100 MeV for the NJL model. In the interior of a star the electron chemical potential is shown 
to be equal to the electric potential energy [46l |47| . In figure [TJd we plot this quantity for stars with 
M = 1.2M for NJL (thick lines) and MIT (thin lines), and for stars with M = 1.4M Q for NJL (thick 
lines). In the last two cases, which correspond to stars close to the maximum mass for a stable star, the 
curves for S = 0, 1 and 2 are almost coincident and so we have only plotted the S = curve. 

One of the motivations of this work was to verify whether, with the NJL model, the presence or absence 
of a crust on a quark star affects its properties and how these properties differ from those obtained with 
the MIT model. We already knew that within the MIT and the CFL models the differences are non- 
negligible [25^. From Table I and Fig.^lthe presence of the crust, described by the BPS EOS affects 
the MIT model than the NJL model. With the NJL model the main differences appear in the radius and 
the central energy density. 

We note that our results are different from those obtained in |48j . where a lower maximum mass was 
obtained, probably due to a different choice of the parametrization for the NJL model. 



7 



TABLE I: Quark star properties for the EOSs described in the text, 
type crust entropy neutrinos M max Mj, max R e 
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IV. CONCLUSIONS 

In this paper we consider the properties of quark stars with the NJL model for the EOS. This model has 
more realistic features than the MIT bag model, in particular, chiral symmetry. Due to the large s-quark 
mass at low densities, several features of the NJL stars are different from the corresponding quantities for 
the MIT stars: smaller strangeness content, higher electron chemical potentials, smaller neutrino fractions 
in stars with trapped neutrinos, higher maximum star masses and smaller maximum central densities. In 
fact, the properties of quark stars with the NJL model are closer to the corresponding properties of hybrid 
stars, e.g. maximum masses, than to the properties of the MIT quark stars. 

We also show that if a quark star is formed directly from a supernova explosion there can be no delayed 
blackholc formation, because the maximum baryonic masses of hot stars are always smaller that the 
corresponding masses of cold stars. 

The electron chemical potential inside the stars has quite large values with the NJL model. This 
influences the properties of the surface of a quark star and the possibility of crust formation. 
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FIG. 2: EOS obtained for S = 0, 1, 2 without neutrinos and with trapped neutrinos within the MIT bag model 
(top figure) and the NJL model (bottom figure). S = gives always the softer and S = 2 the harder EOS. 
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FIG. 4: Strangeness content obtained for S = 0,1,2 without neutrinos and with trapped neutrinos within the 
MIT bag model (top figure) and the NJL model (bottom figure). 
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FIG. 6: Mass-radius plots for bare and crusted MIT and NJL stars at S = with and without neutrinos. 
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FIG. 7: The electron chemical potential a) as a function of density; b) in the interior of stars with M = 1.2 and 
1.4 Mq. Thin lines are for the MIT model and thick lines for the NJL. 



